Dynamic hysteresis in cyclic deformation of crystalline solids 



O 
00 



Lasse Laurson and Mikko J. Alava 
COMP Centre of Excellence, Department of Applied Physics, 
Aalto University, PO Box I4IOO, 00076 Aalto, Espoo, Finland. 

The hysteresis or internal friction in the deformation of crystalline solids stressed cyclically is 
studied from the viewpoint of collective dislocation dynamics. Stress-controlled simulations of a 
dislocation dynamics model at various loading frequencies and amplitudes are performed to study 
the stress - strain rate hysteresis. The hysteresis loop areas exhibit a maximum at a characteristic 
frequency and a power law frequency dependence in the low frequency limit, with the power law 
exponent exhibiting two regimes, corresponding to the jammed and the yielding/moving phases 
of the system, respectively. The first of these phases exhibits non-trivial critical-like viscoelastic 
dynamics, crossing over to intermittent viscoplastic deformation for higher stress amplitudes. 
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The response of interacting many-body systems to os- 
cillating external fields is an old problem in physics, with 
many applications in materials science and engineering. 
In general, due to the competing time scales of the inter- 
nal relaxation and the external perturbation the response 
will generally be out of phase with respect to the exter- 
nal field [Ij]. This gives rise to a dynamic hysteresis loop 
with an area depending on the driving frequency and am- 
plitude. In a magnet driven by an oscillating magnetic 
field h{t), the magnetization m(t) lagging behind the field 
leads to a non- vanishing hysteresis loop area A = f mdh 
0. Hysteresis as such is a very general phenomenon, 
and has been studied in many contexts as from the me- 
chanical response of materials Q, to electronics 0: cell 
biology , neurobiology Q and quantum systems Q . 

Mechanical dissipation or internal friction is one man- 
ifestation of the dynamics of dislocations in crystalline 
solids. Stress-strain hysteresis, in stress or strain con- 
trolled experiments 0], is related via the hysteresis loop 
area to the energy dissipated per cycle. Since disloca- 
tions are line-like objects, internal friction has also been 
described microscopically by the back-and-forth dissipa- 
tive motion of individual dislocation segments [§]. How- 
ever, plastic, irreversible deformation has been shown 
over the last decade to be a highly co-operative process 
with avalanche dynamics and long-range spatio-temporal 
correlations Q. Even the simplest dislocation dynam- 
ics models - which nevertheless describe to a large de- 
gree some real materials - demonstrate phenomena like a 
yielding/jamming transition at an applied stress a — ac 
separating a phase with frozen dislocations from a mov- 
ing phase with a stress-dependent av erage s train rate - 
the order parameter of the transition 'l0l-ll2l|. This is in 
analogy to systems exhibiting criticality due to a depin- 
ning transition - separating in the adiabatic, thermody- 
namic limit frozen and active states, with the order pa- 
rameter given by the average velocity - such as interfaces 
in random media and vortices in type-II supercon- 
ductors [li l. 

In this Letter we consider the dynamic strain rate hys- 



teresis of dislocation assemblies from the viewpoint of 
collective dislocation dynamics. The important aspects 
are i) the various behaviors in the phase diagram, char- 
acterized by the amplitude and the frequency of the ex- 
ternal driving, ii) the collective phenomena that underlie 
the observations from the simulations, and iii) the theo- 
retical and experimental implications of our results. We 
discuss the scaling of the hysteresis, and link it to a pic- 
ture related to depinning transitions. Recent theoreti- 
cal ideas suggest that due to the long-range dislocation 
stress fields this transition should be described by the 
mean field depinning transition [TBI . 16 1. However, our 



results point out that this simple picture is incomplete, 
calling for novel theoretical ideas to properly describe 
the glassy, critical-like dynamics observed in the jammed 
phase of the system. 

Dislocation physics has been recently studied with 
ma ny simp lified models from discrete dislocation dynam- 



17h19{ to phase field [20i| and automaton models [2l| 



Here, we consider the stress-controlled hysteretic dynam- 
ics within a two-dimensional discrete dislocations dynam- 
ics model Such a model captures many of the in- 
teresting aspects of real crystal plasticity, including the 
scale free distribution of avalanches of plastic deforma- 
tion as well as an Andrade primary creep law 
[1^ . It represents a cross section (xy plane) of a single 
crystal with a single slip geometry and straight parallel 
edge dislocations along the z axis. The N dislocations 
glide along directions parallel to their Burgers vectors 
b — ±.hux- Equal numbers of dislocations with positive 
and negative Burgers vectors are assumed, and dislo- 
cation climb is not considered for simplicity. The dis- 
locations interact through their long-range stress fields, 
cr,(f) = Dhx{x'^-y'^)/{x'^ + v'^f, where D = /i/27r(l-i^), 
with jjL the shear modulus and u the Poisson ratio of 
the material. The overdamped equations of motion read 
Xd^Vn/h = s„6E„^„ Smas{rn7n) + <j{t)\, with Vn the ve- 
locity and Sn the sign of the nth dislocation, and Xd is 
the dislocation mobility, implicitly including effects due 
to thermal fluctuations. a{t) is the sinusoidal external 
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FIG. 1. (Color online) Examples of typical "locked-in" strain 
rate hysteresis loops reached by the system after a transient, 
for low and high stress amplitude (left and right column, re- 
spectively) and for various frequencies. All loops exhibit a 
clockwise rotation direction. For small ao and/or large uj the 
loops are smooth, but become increasingly intermittent on 
increasing ao and decreasing cj. Notice that due to the large 
variation of et-values for the different parameter values con- 
sidered, the et axis scales are different in different sub-figures. 

stress, a{t) = tTosin(a;i), with cto the amplitude and uj 
the angular frequency. The equations of motion are inte- 
grated with an adaptive step size fifth order Runge-Kutta 
algorithm, by measuring lengths in units of 6, times in 
units of l/{xdDb), and stresses in units of D, and by 
imposing periodic boundary conditions in both x and 
y directions. Two dislocations of opposite sign with a 
mutual distance smaller than 2b are removed from the 
system (dislocation annihilation). 

The simulations are started from a random initial con- 
figuration of A'o = 1600 dislocations within a square cell 
of linear size L = 2006. The system first relaxes with 
a{t) = 0, to a metastable dislocation arrangement. After 
the annihilations during the relaxation, N = 500 — 600 
dislocations remain. Then, the oscillating external stress 
is turned on, and the evolution of the system is mon- 
itored by measuring the time dependence of the strain 
rate, de{t)/dt = €t{t) = b/L^ J2n ^nVnit)- We simu- 
late the system extensively for a wide range of ctq and 
w-values, with several realizations of the random initial 
configuration considered in each case. 

The resulting stress - strain rate hysteresis loops ex- 
hibit a variety of properties, depending on ctq and tu. 
After an initial transient, the system tends to settle into 
a "locked-in" steady state (usually reached within the 20 
cycles we consider) in which the same hysteresis loop is 



repeated with a clockwise rotation direction in the a-ct 
plane. Fig. [1] shows examples of such locked-in loops 
for different ctq smd uj. For small ctq and large uj (i.e. 
under conditions where a typical distance traveled per 
cycle by a dislocation is small), the loops are smooth 
and the strain rate et{t) obeys sinusoidal dynamics with 
a well-defined phase difference compared to the external 
drive. During the initial transient leading to this smooth 
steady state, the system typically exhibits bursty dis- 
location rearrangements, but will settle into a smooth 
locked-in state after a few cycles. For larger ctq and/or 
smaller w, even the steady state cyclic dislocation dynam- 
ics becomes intermittent, characterized by avalanche-like 
dislocation rearrangements. Interestingly, also in this 
case the system is usually able to find a locked-in steady 
state within the 20 cycles we consider, repeating the same 
bursty dynamics during each cycle in the steady regime. 
The transient time to reach the steady state tends to in- 
crease upon increasing ctq and decreasing w. For large 
ctq and low uj (bottom right corner of Fig. [T]), the loops 
exhibit curvature consistent with the idea that the low- 
frequency limit is described by et ^ (c — ctc)^, with (5 > 1 
0. 

We proceed to characterize the intermittency of the 
steady state cyclic dislocation dynamics, by considering 
the average normalized absolute deviations of et{t) from 
a best-fit sinusoidal function, 

Aet = (|et(t) - et,osin(a;i + a;o)|)/et,o, (1) 

where etfi and are fitting parameters, and (...) in- 
dicates an average over both time and different initial 
configurations. Large values of Aet indicate the presence 
of non-trivial or intermittent dynamics. Fig. [2] shows 
Aet as a function of ctq for various w, demonstrating that 
the intermittency increases with ctq and decreases with 
LJ. By applying a threshold value for Aet, one finds a 
phase boundary separating "phases" with smooth and 
intermittent dynamics in the (Tq — w plane (inset of Fig. 
121) ■ The precise location of this boundary depends on the 
threshold value used, but qualitatively the phase diagram 
looks the same for a range of threshold values. 

Our main result concerns the area Ahyst{<y(),uj) = 
§ etda of the steady state stress - strain rate hystere- 
sis loops as a function of cto and uj. These are summa- 
rized in Fig. [31 Ahyat{o'Q^uj) exhibits a maximum at a 
characteristic frequency uj* w 0.06 independent of (Tq, 
corresponding to the resonance frequency of the effective 
confining potential (see the oscillator model below). The 
A/iysf (a;)-data for various uo can be collapsed by nor- 
malizing with (Tg, leading to two distinct low frequency 
power laws Ahyst ~ w'^, with exponents k « 0.82 and 
K « 0.69 for do < 0.015 and ctq > 0.015, respectively. 
The stress amplitude value ctq = crc{N) « 0.015 separat- 
ing these two regimes corresponds roughly to the maxi- 
mum (T-value for which the power-law Andrade creep is 
observed in a constant stress simulation, i.e. et ^ 
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FIG. 2. (Color online) Normalized absolute deviations of et(t) 
from a best-fit sinusoidal function, et,o sm{ijjt + iJo), character- 
izing the intermittency of the cyclic dislocation dynamics in 
the "locked- in" steady state as a function of ctq, for various uu. 
The inset shows a phase diagram in the ui-ao space, displaying 
smooth and intermittent phases separated by a phase bound- 
ary obtained by thresholding the data in the main figure, with 
the threshold shown as a dashed line. 



with 8 close to 2/3 [H [H [li-ii]. For a larger applied 
stress in the DC-driven case, the system would reach a 
(quasi)stationary moving/ flowing state with a non-zero 
mean strain rate [lo| . Thus, we argue that the two stress 
amplitude regimes with the different K-values correspond 
for the system to in the jammed {ctq < adN)) and mov- 
ing (ctq > ac{N)) states for a constant external stress. 

We note that the magnitude of Ahyst is related to both 
the phase difference between cr(t) and et{t), and to the 
strain rate amplitude, i.e. loq and et^o in Eq. ([1]). Fig. 
|4] shows that for (Tq < cTc, wq is independent of ctq and 
approaches 7r/2 for w — ?> 0, and goes to zero for large 
(jj. Viscoelasticity is typically characterized by the phase 
lag S between a{t) — ctq sin(a;i -I- S) and e{t) = eo sin(a;i), 
with 6 = and 7r/2 corresponding to perfectly elastic 
and viscous dynamics, respectively. Thus, the relation 
'j-'o = 7r/2 — (5 implies that the dynamics extrapolates 
between perfect elasticity for a; — > and perfectly viscous 
dynamics in the high-frequency limit. For (Tq > Uc, ujq 
starts to decrease for small w, indicating the presence 
of plastic dislocation rearrangements, also visible in the 
intermittency of the dynamics (Fig. [S])- Rescaling the 
strain rate amplitude by ctq leads to a data collapse for 
Co < CTc, with et.o/co w'', k « 0.82 for small cj, while 
for large uj, et,o/<^o Nb/L"^ « 0.0125, corresponding 
to N ^ 500 dislocations freely following a{t) in a system 
of size L = 2005. For (Tq > (Tc, there are deviations 
from the low-frequency power law, again corresponding 
to intermittent viscoplastic deformation. 

As a naive attempt to understand the observed scaling 
behavior, we consider a simple oscillator model for the 
dislocation dynamics. In general, dislocations will oscil- 



FIG. 3. (Color online) The scaled hysteresis loop areas 
Ahyst/ o-Q for various cro as a function of ijj, showing a data col- 
lapse with two distinct low-frequency power laws, Ahyst ~ tij", 
with K ^ 0.82 for cro < 0.015 (Tc and k. f» 0.69 for 
(Jo > 0.015. The inset shows the measured k as a function of 
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late back and forth due to the sinusoidal applied stress. 
However, dislocation interactions induce a tendency to 
form dislocation structures of varying complexity - dislo- 
cation multipoles - with each multipole moving together 
with a strain rate in a way dictated by its net Burgers 
vector under the applied field and interactions with 
the rest of the system. For small ctq and large w we de- 
scribe the latter by a harmonic potential, and write the 
equation of motion for the ith multipole as 

^L^Ml. ^ it) = -K(^) [x(') {t) - x(') (0)] -t- ao sin(wi), 

where K^"^^ is the effective spring constant characterizing 
the confining potential of the ith multipole (with a center 
of mass x^'^\t)) due to long-range interactions with the 
other multipoles. The asymptotic {t oo) solution of 
Eq. ^ is given by 



(i) fo^'Vo [K^^'^io cos(wt) -I- sin(wt)] 



(3) 



The total strain rate is obtained by summing over the 
multipoles, et = X^i^t'^- Disregarding fiuctuations by 
setting = beff and K'-'^ = Keff for all i, one obtains 

beffNmp<Jo KeffUJCOs{u!t) +uj'^sm{ujt) 
- 7^. -fPi — r-3. ' (4) 



L2 



with Nmp the number of dislocation multipoles in the 
system. Notice that Eq. (|4|) corresponds to a clockwise 
direction of rotation in the a — tt plane, as observed in 
the simulations. The area of the hysteresis loop, Ahyst = 



4 



etda, is given by 



Ahystil^, CTo) = 



beffNmpfTQ nKeffU} 



^eff 



(5) 



([5]) predicts a maximum of A^yst around uj = 
i Keff, and a power law frequency dependence 



and Ahyst ^ 
Fitting Eq. 



for oj ^ UJ* and w -C w*, 
© to the data in Fig. [3] 



Eq. 

Afiyst U! ^ 

respectively. 

leads to beffN^p ~ 400 independent of (Tq, suggesting 
that most dislocations would move either individually or 
within wall-like structures. However, while such a sim- 
ple model results in the stress amplitude dependence ob- 
served in simulations, i.e. Ahystio'o) ~ Cg, it obviously 
fails to reproduce correctly the non-trivial low-frequency 
K-exponents. 

Thus it is necessary to go beyond such simplistic de- 
scriptions by considering ideas from critical phenomena, 
applied to a yielding transition 12 1. It has been 
proposed that due to the long-range dislocation stress 
fields, this should be described by the mean- field depin- 
ning transition 15|, |16|. To test this idea within the 



present framework, we proceed to contrast our results 
with those obtained recently for a mean field elastic in- 
terface subject to AC driving [2^. There, the exponent 
of the low- frequency power-law Ahystito) ~ w'^ describing 
the force-velocity (the latter being the order parameter 
of the depinning transition) hysteresis loop area has been 
shown to exhibit three regimes as a function of the ap- 
plied force amplitude (Tq, with the exponent k assuming 
the values k ~ 0.67 or 0.75 (for cusped and smooth disor- 
der, respectively f26|) for ao ^ ac, n ~ 0.82 for ctq ~ Cc 
and K « 1 for (Tq ^ ctc- <^c is the critical depinning 
force of the DC-driven system. Thus, in particular, our 
numerical results do not agree with the mean field depin- 
ning scaling of the loop area for small force / stress ampli- 
tudes, corresponding to the pinned/jammed phase: The 
pinned phase of the mean field interface exhibits trivial 
dynamics, with the K-exponent coinciding with that of 
the naive oscillator model, whereas we find here a k < 1 
for a wide range of stress amplitudes with ao < ac- An 
additional difference is that for dislocations we observe 
only a single hysteresis loop [27|, while for interface de- 
pinning models one typically observes a secondary loop 
with counterclockwise rotation direction for (Tq ^ ac in 
the region cr > 0-^ [ll 111 . 

Consequently, our results reveal that various scaling 
features in the dynamic hysteresis of crystalline solids 
relate to the collective dynamics of dislocations. From 
the theoretical point of view the central observation is 
that unlike the pinned phase of conventional AC-driven 
mean field interfaces, the jammed phase of the disloca- 
tion system exhibits critical-like dynamics. In fact, the 
K-value we find for ctq < cFc coincides with the mean 
field depinning result for ctq ss ctc (k ~ 0.82), suggesting 
that the system exhibits criticality in the entire region 
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FIG. 4. (Color online) The phase difference uiq (upper curves) 
and the scaled strain rate amplitude et,o/o"o (lower curves) 
for various ao. For ao < (Jc, cjo approaches 7r/2 for tj — > 
(dashed line, a signature of perfectly elastic dynamics), while 
for larger ao, cjo starts to decrease for small ui, indicating 
the presence of plastic rearrangements. et,o/ao obeys a power 
law et,o/ao ~ 1^" for ao < ctc (solid line, with deviations 
from the power law for ao > ac), and approaches a value 
et,o/ao = Nh/L^ ~ 0.0125 (dash-dotted line) for large u. 



Q < < Uc- Similar observations have been made by 
Ispanovity et al. 



25|, who found that the dislocation 
system subject to a small constant stress exhibits glassy 
power-law relaxation up to a time scale limited only by 
the system size rather than the applied stress value. We 
think this is due to the dynamic nature of the effective 
disorder - the dislocations are subject and jam due to 
a random stress field rather than to quenched disorder 
absent here, but fundamental to conventional depinning 
models. Interesting extensions of the present study could 
include considering the effect of a non-linear mobility law 
P, '^i', ^2^ (which might give rise to a dynamic transition 
liiLSOj), or cyclic loading of alloys exhibiting the Portevin- 



31|,|32| 



Le Chatelier effect 

To conclude, the time-dependent loading of disloca- 
tions exhibits features that result from collective dynam- 
ics. This suggests that such signatures should be seen 
during the deformation of any material containing dislo- 
cations, and that they should also be looked-for in the 
yielding of non-crystalline materials 33|. These findings 
call for new experimental (for instance in colloidal crys- 
tals ^]) and numerical studies (e.g. molecular dynamics 
simulations [s^l) of cyclic dislocation dynamics, as well 
as novel theoretical ideas to properly describe dislocation 
jamming. 
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